Particle Scale Dynamics in Granular Impact 



Abram H. Clark, 1 Lou Kondic, 2 and Robert P. Behringer 1 

department of Physics & Center for Nonlinear and Complex Systems, Duke University, Durham, NC 27708, USA 
Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, NJ 07102, USA 

(Dated: October 26, 2012) 

We perform an experimental study of granular impact, where intruders strike 2D beds of photoe- 
lastic disks from above. High-speed video captures the intruder dynamics and the local granular 
force response, allowing investigation of grain-scale mechanisms in this process. We observe rich 
acoustic behavior at the leading edge of the intruder, strongly fluctuating in space and time, and we 
show that this acoustic activity controls the intruder deceleration, including large force fluctuations 
at short time scales. The average intruder dynamics match previous studies using empirical force 
laws, suggesting a new microscopic picture, where acoustic energy is carried away and dissipated. 
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The penetration of a dense granular material by a high- 
speed intruder occurs routinely in meteor and ballistic 
impacts. Many previous studies jlJ-Q , both recent and 
dating back to Euler and Poncelet, have used variations 
of a macroscopic force law: 

F = mz = mg — f(z) — h(z)z 2 . (1) 

Here, z is the intruder depth relative to the top of the 
original, unperturbed surface (i.e., z = at initial im- 
pact), mg is the of gravity force, f(z) characterizes hydro- 
static effects, h(z) is often assumed constant, h(z) = b, 
and dots denote time derivatives. In Eq. (p}, h(z)z 2 rep- 
resents a coarse grained collisional stress. We note that 
other effects, including a depth-dependent Coulomb fric- 
tion term have been proposed la, [6f| . Despite the success 
of extensive previous studies [iHla] , the connections be- 
tween the local granular response, the microscopic pro- 
cesses responsible for dissipating kinetic energy, and the 
dynamics of the intruder are still subjects of debate, 
largely due to experimental difficulties in obtaining suf- 
ficiently fast data at small scales. 

In this Letter, we address this issue experimentally by 
high-speed imaging of an intruder of mass, to, which 
impacts a quasi-two-dimensional system of photoelas- 
tic particles (bidisperse, larger particle diameter d) at 
speeds vq < 6.5 m/s, yielding both the intruder dynam- 
ics and the force response of individual grains (Fig. Q]). 
Here, as in many previous experiments, v << C, where 
G ~ 300 m/s is the granular sound speed, measured 
from photoelastic space-time plots, as in Fig. QJb). The 
frame rates of ~ C/d capture the microscopic granular 
response. The primary intruder energy loss mechanism in 
these experiments is due to intense, intermittent acoustic 
pulses traveling at speeds ~ G along networks of grains, 
transmitting energy from the intruder into the medium. 
These pulses decay roughly exponentially with distance 
from the intruder. The force on the intruder is strongly 
fluctuating, due to the intermittency of the force net- 
work/acoustic activity, but the mean behavior is consis- 



tent with empirical models used previously [lHS|. 

Experimental Techniques. -The experimental appara- 
tus consists of two thick Plcxiglas sheets (0.91 m x 1.22 m 
x 1.25 cm), separated by a thin gap (3.3 mm) which is 
filled by photoelastic disks (thickness of 3 mm) of two 
different diameters (6 mm and 4.3 mm). These disks are 
cut from PS-1 material, (Vishay Precision Group; bulk 
density of 1.28 g/cm 3 , elastic modulus of 2.5 GPa, Pois- 
son's ratio of 0.38). Intruders are machined from bronze 
sheet (bulk density of 8.91 g/cm 3 , thickness of 0.23 cm) 
into disks of diameters D of 6.35 cm, 10.16 cm, 12.7 cm , 
and 20.32 cm (data for D — 12.7 cm intruder used in im- 
ages and time-series data shown here is typical for all D). 
These intruders are dropped from a height, H < 2.2 m, 
through a shaft connected to the top of the thin gap 
containing the particles, producing an impact speed, 
vo ~ {2gH) 1 / 2 . Results are recorded with a Photron 
FASTCAM SA5, at a resolution of 256x584 pixels (~ 10 
pixels per d), at 40,000 frames per second. To locate and 
track the intruder, we use a circular Hough transform at 
each frame. Velocity, v, and acceleration, a are calcu- 
lated by numerical differentiation, with a low-pass filter, 
cut-off frequency of 133 Hz ~ (7.5 ms) _1 ~ v /D), ap- 
plied with each derivative to reduce noise amplification. 
The frequency cutoff is as large as possible while main- 
taining a signal-to-noise ratio of 10:1. This yields inter- 
mediate time scale data for v (vi nt ) and a (dint) which 
are still strongly fluctuating in time. Photoelastic im- 
ages are normalized by a calibration image, taken before 
the intruder is dropped, to account for imhomogeneties in 
the light source. After this, the discrete gradient-squared 
(G 2 = | V/| 2 ) of the image is computed using the spatial 
variation of the image intensity, J; the sum of the G 2 in 
a particular region measures the local force response [3] 
(i.e. beneath the intruder, as in Fig. [TJ). A static cali- 
bration, covering the full range of G 2 encountered in any 
impact, was performed by placing a weighted piston on 
a box of about 100 particles that are subject to the same 
light intensity as used in the experiments. As shown in 
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FIG. 1: (color online), (a) Six selected frames, starting at 2.75 ms after impact and spanning 475 /is, showing the end of a 
typical compressional event which generates an acoustic pulse, which disconnects from the intruder, (b) A space-time plot of 
G 2 in an angular region under the bottom half of the intruder (half-annulus) over time. The x-axis is time, and the y-axis is 
radial distance from the bottom of the intruder, where the top of the plot corresponds to the bottom edge of the intruder. The 
slope of the disturbances gives a consistent acoustic speed of ~ 325 m/s. (c) The sum of the response in the space-time plot 
above, after subtracting background inhomogeneities. Calibrating this will yield our measurement of instantaneous force, as 
shown later, where the range shown above (0 to 1.1 arbitrary units (AU)) maps to an intruder acceleration range of to 27 g. 



Fig. [2fb), it is essentially linear. 

Comparing to previous models. -An important question 
is whether the observed dynamics are consistent with ex- 
isting models, i.e. Eq. JT]). To address this, we consider 
the intruder trajectory, z(t), and the filtered derivatives 
Vi n t and ciint- As noted, the derivatives, particularly 
ciint, are strongly fluctuating, and these fluctuations are 
a physical aspect of the dynamics, as discussed below. 
Plots of Oint vs. v 2 nt data from different impacts with 
varying vq, show good agreement, within fluctuations, 
with Eq. {!]). This analysis yields f(z) and h(z): a con- 
stant value for h(z) (i.e. h(z) = b ~ 5D) after an initial 
transient at impact, and f(z) which is nearly linearly 
increasing in depth. 

However, for any individual trajectory, we measure 
large fluctuations in a, n j (Fig. [2]), on a scale that is com- 
parable to the mean acceleration. These fluctuations are 
absent in the "slow-time" models discussed above, and 
their large amplitude is both a novel observation and a 
potential weak point of the models. That is, the braking 
of the intruder is not a smooth steady process, but a se- 
ries of events where the intruder is subjected briefly to 
large accelerations, followed by more quiescent periods 
that can be close to acceleration-free. 

Connecting acoustic activity to intruder deceleration. - 
As noted, during an impact, we observe complex propa- 
gating force networks (known as force chains) generated 
intermittently at the leading edge of the intruder as it 
moves through the medium, as shown in Fig. HJa), as 



well as in Supplementary Videos 1 and 2. To quantify 
the photoelastic response, we consider the angular re- 
gion extending radially outward from the bottom half of 
the intruder over a length ~ 10c?, forming a half-annulus. 
Fig. [lib) shows a space-time plot of the total photoelastic 
response in this region. 

To relate the photoelastic activity to the acceleration 
fluctuations, we compare the total photoelastic response 
in the angular region immediately under the intruder to 
dint (Fig. [2]). Photoelastic data is obtained at 40 kHz, 
which is about 500 times faster than the frequency cut- 
off for ai n t. Comparing ai n t to the photoelastic response, 
G 2 , requires time-filtering the photoelastic data such that 
the time scale matches that of eti„ t . This gives a com- 
parison at the intermediate time scale; a plot, Fig. [2 of 
Oint and filtered G 2 data gives the same curve, showing 
that the two are virtually identical. For this comparison, 
we first normalized G 2 by a constant to obtain the op- 
timum agreement between filtered G 2 and ai n t, but this 
normalization matches well with the static calibration of 
G 2 discussed above. (We used this double comparison to 
be sure that the static calibration matched well with the 
dynamics measurements.) We conclude that the large 
photoelastic events are the main force mechanism acting 
on the intruder. By inference, the energy loss for the in- 
truder is tied to these acoustic events, rather than, e.g., 
to frictional drag with the intruder. 

Acoustic dissipation. -Once the acoustic pulses have 
moved ahead of the intruder, there must be a loss mecha- 
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FIG. 2: (color online.) (a) Comparing the intruder trajec- 
tory to photoelastic response from Fig. \V[c) shows that the 
intruder acceleration is very well correlated to the photoe- 
lastic/ acoustic fluctuations in high-speed videos. We time- 
average the photoelastic response (thick, blue line) to match 
the time scale of the acceleration measurement, dint (black, 
dashed line), which has limited time resolution. Rescaling 
the photoelastic measurement gives extremely close agree- 
ment with the measured deceleration (both the mean and 
fluctuations). The calibrated photoelastic force measurement 
without time-filtering (thin, red line) shows much larger fluc- 
tuations at a much shorter time scale, (b) Inset shows calibra- 
tion of photoelastic response vs. 2D pressure (force per width 
of intruder or piston) from experiment (black dashed line) and 
from a static test (blue circles), with good agreement. 



nism of these disturbances within the material. Hence, it 
is important to examine how fast and how far the acoustic 
pulses propagate. To this end, we observe the photoelas- 
tic response in a long, thin angular slice, centered directly 
beneath the intruder with a width of 7r/8, which extends 
25d beneath the intruder. Space-time plots of the re- 
sponse in this region indicate a wave speed of about 325 
m/s (~1/10 of the sound speed in the bulk material from 
which the particles are cut). To determine the attenua- 
tion of the acoustic pulses, we plot normalized intensity 
versus depth. The normalization for each pulse is the cu- 
mulative photoelastic response, G 2 , over its full duration. 
The normalized photoelastic response averaged over mul- 
tiple events shows an exponential decay (Fig. [3]), with a 
decay length of ~10 particle diameters, which is short 
enough that reflections from the bottom or sides of the 
container are not important. It is unclear which grain- 
scale interactions are responsible for this decay, but it 
could be explained by force-chain splitting, grain-grain 
friction, restitutional losses for each "collision" , or other 
dissipative mechanisms. 

Fluctuation statistics and stochastic description. - 
Large fluctuations in the photoelastic response (Fig. [3]) 
suggest a stochastic description, which captures mean be- 
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FIG. 3: (color online.) Photoelastic pulses decay as they 
propagate away from the intruder. We observe a thin an- 
gular slice of opening angle of tt/8 radians, extending 25d 
below, and centered directly beneath the intruder. We use 
40 different pulses from different impacts of a single intruder 
(D — 21.17d), where the intruder velocity at the pulse emis- 
sion varies between 2 and 6 m/s. We then plot the natural 
logarithm of G 2 per area as a function of depth for each pulse, 
normalized by the total intensity in the pulse (wave intensity 
will decrease as 1/r moving away from a point source in 2D, 
and this effect has already been accounted for in this plot). 
Imposed fit (thick, red line) is exp(— r/L), where L is the 
decay length, roughly 10 particle diameters. 

havior as well as short-time fluctuations. For example, 
one might modify Eq. (TJ} to: 

F(z, z, t) = mg - [f(z) + h(z)z 2 } V (t). (2) 

Here, 77 (t) is a multiplicative stochastic term, which 
should follow directly from microscopic physics and have 
a mean of unity. A multiplicative term is chosen here 
since rescaling by the mean photoelastic behavior yields 
a statistically stationary fluctuating term, as discussed 
below, and since fluctuations in dense granular systems 
often scale with the mean (as here) . 

To experimentally characterize the fluctuations in 
Eq. ©, we write 77(f) - G 2 (t)/G 2 avg {t), where G 2 (t) is 
the photoelastic time series used to measure force (e.g. 
bottom of Fig. [1]), and G 2 vg (t) is the mean behavior, ob- 
tained by fitting a low-order polynomial to G 2 (t). This 
yields a fluctuating term which appears statistically sta- 
tionary throughout the duration of an impact, as shown 
in Fig. |U Typically, r)(t) has an autocorrelation decay 
time of ~1 ms, and a probability distribution function 
(PDF) that is nearly exponential. The PDF describes 
the likelihood of the large events which dominate the de- 
celerating force. Such a PDF is typical for forces in static 
dense granular systems and is presumably related to the 
probability of generating force-chain-like structures. 

Surprisingly, the fluctuation statistics show almost no 
dependence on intruder size. One might expect that the 
contact forces or force chains generated from two suf- 
ficiently separated points along the bottom of the in- 
truder are uncorrelated. If so, increasing the intruder size 
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FIG. 4: (color online.) (a) The fluctuating term, r)(t), 
for a single impact (typical for all impacts), where n(t) ~ 
G 2 (t)/G 2 avg {t), as discussed in the text, (b) The autocorrela- 
tion and (c) the PDF of the combined fluctuating signals for 
all impacts for each intruder (~20 runs per intruder). The 
semi-log PDF plot shows P{rj) ~ exp(— 77). We see an auto- 
correlation time of ~ 1 ms, which gives a typical event time, 
which agrees with video frames in Fig. Ufa). Lack of depen- 
dence on intruder size suggests a collective mechanism and not 
a simple combination of uncorrelated, random force chains. 



would include more of these independent forces, which, 
by the central limit theorem, would yield smaller and 
more Gaussian-like fluctuations, regardless of the statis- 
tics of each one. However, this does not occur, suggesting 
a more subtle collective mechanism. One possibility is 
that spatially separated intruder-particle contacts often 
excite the same persistent force network. 

Conclusion. -hi this Letter, we present a new micro- 
scopic picture of the force on an intruder moving through 
a granular material, which focuses on acoustic activity 
and fluctuations due to the generation of force-chain-like 
pulses. We observe consistency with established impact 
force models, but with substantial fluctuations in the 
measured deceleration of the intruder during the impact 
process. We have shown that the acceleration profiles, 
including these fluctuations, are a direct consequence of 
acoustic pulses transmitted along networks of particles. 
Other recent studies have indicated an important role 
for granular force networks in intruder impacts 17| and 
acoustic transmission The microscopic description 

presented here should also help connect granular impact 
experiments with differing microstructure, such as more 
dilute or compacted [8|, [l9j, or anisotropic (e.g. sheared) 



systems, or even more general experiments on granular 
flow around an obstacle. Strong force fluctuations sug- 
gest a stochastic model, which gives a natural way to sep- 
arate the slowly varying macroscopic response from fast- 
time fluctuations. We believe that the granular sound 
speed is critical in our description, so we expect substan- 
tial differences when intruder speeds arc close to sonic 
or even supersonic. This could be achieved by increasing 
intruder velocity or reducing the granular sound speed by 
using softer material. Also of interest is how these effects 
translate to three-dimensional systems, or systems with 
much larger ratio of intruder size to particle size. 

This work has been supported by the US DTRA un- 
der grant HDTRAl-10-0021. We very much appreciate 
additional input from Drs. C. O'Hcrn and W. Losert. 
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